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1. INTRODUCTION 

The temperature and moisture regimes in a forest are key components in the forest ecosystem dynamics. 
Observations and studies indicate that the internal temperature distribution and moisture content of the tree 
influence not only growth and development, but onset and cessation of cambial activity[l], resistance to insect 
predation[2], and even affect the population dynamics of the insects[3]. Moreover, temperature directly affects the 
uptake and metabolism of pollutants from the soil into the tree tissue [4]. Additional studies show that soil and 
atmospheric temperatures are significant parameters that limit the growth of trees and impose treeline elevation 
limitation[5]. 

Directional thermal infrared radiance effects have long been observed in natural backgrounds [6]. In earlier 
work, we illustrated the use of physically-based models to simulate directional effects in thermal imaging[7-8]. In 
this paper we describe the simulated spatial and temporal patterns of soil moisture, temperature, and net radiation 
in a deciduous forest using our recently developed three-dimensional, macro-scale computational tool that 
simulates the heat and mass transfer interaction in a soil-root-stem system (SRSS). The SRSS model includes the 
coupling of existing heat and mass transport tools to simulate the diurnal internal and external temperatures, 
internal fluid flow and moisture distribution, and heat flow in the system. 

2. METHODS 

We simulate a tree and its root system in the SRSS as three different functional materials; the bark, xylem 
(sapwood), and heartwood. The bark material is the outer covering of the tree consisting of dead phloem cells, 
which are air-filled sieve-tube cells that are thermally less conductive than the other materials in the trunk or stem 



and hydraulically non-conductive. The second material, the xylem, is the living, hydraulically active part of the 
wood that is responsible for the transport of fluid (water and nutrients) from the roots into the rest of the plant. 
The innermost material is the heartwood. The heartwood provides structural support to the tree and is composed 
of dense dead xylem cells that are typically not hydraulically conductive. The root system is derived from 
measurements of excavated roots of a tree system (Figure la). The tree and root system is embedded in a soil that 
has properties of typical sandy silt (Figure lb). A diumally varying flow function is imposed on the upper 
boundary of the stem that was derived from an average of fluid velocities measured in-situ. The three 
dimensional distribution of water and temperature in the soil is accomplished by approximating Richards’ 
equation[9] and heat conduction and convection equations, respectively. Soil properties include porosity, 
hydraulic conductivity, pressure heat, thermal conductivity, spectral emissivity and reflectivity. Surface heat 
exchange of the soil includes net radiative energy, sensible heat, latent heat, and heat from precipitation events. 

Radiative heat transfers are calculated by the Monte Carlo method. Energy flux is emitted, scattered, and 
reflected around the system until the energy is either absorbed or exits the system. At the short wavelengths, the 
source term is primarily solar illumination, while at the longer wavelengths; the energy is emitted by all surfaces 
in the system. Two modes of the Monte Carlo method are used here, forward Monte Carlo method for longwave 
radiation and backward Monte Carlo method for shortwave radiation. The two methods are applied in this 
manner primarily for computational benefits. In the case where there are a large number of emitting surfaces, the 
forward Monte Carlo method is the most computationally efficient and conversely in the case of a few emitting 
surfaces, the backward Monte Carlo method is more efficient[10]. 




Fig. 1 . Photograph of the excavated root system and corresponding computational domain 



3. RESULTS 


The simulation was executed on a Cray XT4 computer system using 1072 processors and was driven with 
measured meteorological data collected 2-3 May 2010. The simulation was setup with a measured soil 
temperature depth profile and the soil saturated with fluid by specifying the water table to be just below the 
surface of the soil. Figure 2 illustrates the cross-section sap velocity, distribution of surface net radiation, and 
unsaturated soil region. Figure 3 illustrates the resulting growth of the unsaturated soil region in the soil caused 
by the transpirational uptake of the tree and root system. The desaturated region begins at the surface of the soil 
and gradually increases across the surface and down toward the root system. Flow velocities remain constant 
until the unsaturated region surrounds the root. At that point the hydraulic conductivity of the soil decreases and 
the flow rate decreases for that part of the root structure but at the same time flow rates increase for roots still in 
saturated soil regions. Satellite observations offer valuable tools for land remote sensing of global change [11-12], 
but the highly dynamic and spatially varying nature of the radiant fluxes illustrated here present challenges in 
appropriately interpreting these observations. 



Fig. 2. Combined results of the spatially varying soil moisture, temperature, fluid velocity, and surface net 

radiation. 
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Fig. 3. Growth of the unsaturated soil region due to fluid uptake by the root system from 1000 hrs to 

1200 hrs. 
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